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494 SUPPLEMENTARY MATERIAL 

495 Model construction 

496 Since it is well known that birds see a wider colour spectrum than humans (Cuthill, 2006), we 

497 calibrated the colours printed on the models so that they stimulate the photoreceptors of 

498 potential bird predators in a manner comparable to the real wing colours. The method followed 

499 is similar to that of Cuthill et al. (2005) and Stevens et al. (2007). In brief, the reflectance spectra 

500 of each colour patch on real wings was measured with an Ocean Optics USB4000 spectrometer 

501 (Dunedin, FL), with illumination by a PX-2 pulsed Xenon lamp with a narrow-ended (1/8") probe, 

502 held at a constant distance and a 45 degree angle to the butterfly wings. Three replicate 

503 measurements were taken per colour patch, in each of four specimens selected per colour 

504 pattern. Patches of the same colour on fore- and hindwings were measured separately. Since the 

505 reflectance spectra measurements were highly consistent within each colour patch, between 

506 patches of the same colour, and even between individuals with either the same or different 

507 colour pattern, the averages of all measurements per colour were used for further analysis 

508 (Figure SI). The next step was to calculate the photon catch values (as Endler & Mielke, 2005) 

509 for each of the four single cone types (used in colour vision) and the double cones (used in 

510 achromatic vision) for a typical insectivorous bird (a blue tit, Cyanistes caeruleus), using data 

511 from Hart et al. (2000) and a standard D65 irradiance spectrum. Although neither irradiance 

512 spectra from the field locations nor sensitivity data for the real species of bird predators were 
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513 available, the level of error associated with those factors should be relatively small as different 

514 irradiance spectra have little effect on the outputs of most visual models in terrestrial systems 

515 (e.g. Langmore et ah, 2009) and birds seem to have relatively constrained visual sensitivity 

516 across species (Cuthill, 2006). The colours to be printed in the final models were calculated by 

517 printing a range of colour values around those observed on real wings and measuring the 

5 1 8 photon catch values produced. An iterative process was used, whereby this procedure was 

519 repeated, changing the closest colours slightly until photon catches were as close as possible to 

520 those of real wings (see main text, Figure 3). Final colours were selected based on the criteria 

521 that when printed they produced photon catch values for each cone type that fell within the 

522 range of values for the corresponding colour patches on the real butterflies (as with past work 

523 on camouflage: e.g. Cuthill etal, 2005, and mimicry: Merrill etal. 2012). The calibrated colours 

524 were used to replace the colours on digital images of real wings (see main text), using Adobe 

525 Photoshop®, and these final print outs were used to build the models for the experiment. 

526 Control of potentially confounding factors 

527 Table SI lists all potentially confounding factors and the strategies applied to control for their 

528 potential effects during the experiment to test for frequency-dependent selection. 
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529 Table SI. Potentially confounding factors and applied controls. 



Potentially confounding factor 


Applied control 


Observer bias: When marks were weak some 


Sets of models were assigned randomly to 


observers could consider a model as marked 


observers, making sure that no single observer 


whereas others as intact. 


checked the same set of models twice in a row. 




This way, any observer bias is expected to 




equally affect all observations across all patterns 




and all zones. 


Seasonality: The experiment lasted for three 


Experimental trials were alternated between 


months (July-September 2010). During this 


zones. For instance, the first trial was performed 


period, seasonal variations in rainfall and/or 


in one of the parental zones, the second trial in 


temperature could affect predation patterns. 


the hybrid zone, the third one in the other 




parental zone, and the fourth one back in the 




first zone. Any possible effect of seasonality in a 




particular zone is expected to equally affect all 




tested colour patterns. 


Microhabitat heterogeneity: Predation patterns 


Models of each tested colour pattern were placed 


could vary depending on microhabitat 


one after the other along the sampled transects. 


characteristics such as luminosity or vegetation 


This way it was ensured that all three patterns 


coverage. For instance, predation inside a 


were equally exposed to all habitats where the 


covered forest is likely to differ from that along a 


experiment was carried out. Any effect due to 


forest edge or a river shore. 


microhabitat differences is expected to affect all 




colour patterns equally. 
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new paper wings for each experimental trial. 


paper wings was used for all colour patterns in 


Therefore, paper wings were reused up to three 


all studied zones. Any affect of wings wearing is 


times. Wearing of paper wings could affect the 


expected to equally affect all colour patterns in 


probability of the models being attacked. 


all zones. 
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530 Figure SI 
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SUPPLEMENTARY FIGURE LEGEND 



533 Figure SI (A) Average reflectance spectra of colour patches on real wings in the studied 

534 phenotypes of Heliconius erato. Notice the UV reflectance peak on the white patches of the 

535 highlands race, and to a lesser extent on the yellow patches of the hybrid and lowlands race. (B) 

536 Photos of wings illustrating the UV reflectance on the white patches of the highlands race. 
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